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£ ■ 1 Introduction 



We use the term "integrability" in the meaning that a system or equation under consideration 
possesses a Lax representation or a zero curvature representation. Such systems can be solved 
by the inverse spectral transform method (1ST) [H [2] . Exactly integrable evolution systems are 
£T) ' of interest both for mathematics and applications. In particular, systems of the following form 

in ' 

ut = Uxxx + F(u,v,u x ,v x ,u xx ,v xx ), v t = av xxx + G(u,v,u x ,v x ,u xx ,v xx ), (1.1) 



CN ' where a is a constant, excite great interest since about 1980. The paper [3] is devoted to construc- 

tion of systems of the form (jl.lj) among others. Nine integrable systems of the form (jl.ip and 
their Lax representations have been obtained in the paper. In particular, it contains a complete 
list of three integrable systems (jl.ip satisfying the conditions a(a — 1)^0 and ord(F, G) < 2. 
Here ord = order, ord f < n means that / does not depend on u n , v n , u n+ i,v n+ i, .... Here and 
in what follows, the notations u n = d n u/dx n , v n = d n v/dx n are used. 



Two of the three mentioned systems can be written in the following form 

Ut = U3 + V Ml , v t = — \vz + UU\ — VVl, (1.2) 

u t = u-s + vui, v t = -\vz - vv\ + u\ (1.3) 

and the third system is presented below (see (|3.46a|) ). System (|1.2|) was found independently 
in [3] and the soliton solutions were constructed there. This system is called as the Drinfeld- 
Sokolov-Hirota-Satsuma system. 

This paper contains two results: (i) a list of integrable systems of the form (jl.ip with smooth 
functions F, G and a = —1/2; (ii) differential substitutions that allow to connect any equation 
from the list with JT2J) or (fl~3]) . 

There are many articles dealing with integrable systems, but some of them (see, e.g., (5j El [7]) 
consider multi-component systems. Other papers (see, e.g., [U [U [TUl [HI [12]) contain two- 
component systems reducible to a triangular form. The triangular form is briefly considered 
below. There was possibly only one serious attempt [T3] to classify integrable systems of the 
form (jl.ip using the Painleve test. Unfortunately, fifteen systems presented in [13] contain a large 
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number of constants some of which can be removed by scaling and linear transformations. Note 
that there are two-field integrable evolutionary systems Ut = Au-j + H(u, U\, 112) with a non- 
diagonal main matrix A. For example, an integrable evolutionary system with the Jordan main 
matrix is found in [14j . 

Moreover, about 50 two-field integrable systems of the form (jl.ip with a = 1 can be extracted 
from papers [El [I6j [T71 [TBI [19] that deal with vector evolutionary equations. 

Partial solutions of the classification problem for a = and ord G ^ 1 have been obtained 
in [20] , and in [21 j for divergent systems with a ^ 1. A complete list of integrable systems of the 
form (jl.ip does not exist today because the problem is too cumbersome and the set of integrable 
systems is very large. 

Our tool is the symmetry method presented in many papers. We shall point out pioneer or 
review papers only. In [22] the notions of formal symmetry and canonical conserved density for 
a scalar evolution equation are introduced. These tools were applied to classification of the KdV- 
type equations in [23j . A complete theory of formal symmetries and formal conservation laws for 
scalar equations has been presented in [24]. A generalized theory was developed for evolutionary 
systems in [25]. Review paper |26] contains both general theorems of the symmetry method 
and classification results on integrable equations: the third and fifth order scalar equations, 
Schrodinger-type systems, Burgers-type equations and systems. Review paper [27] is devoted to 
higher symmetries, exact integrability and related problems. Peculiarities of systems (jl.ip have 
been discussed in [21] . For the sake of completeness, the main points of the symmetry method 
and some results necessary for understanding of this paper are considered in the Sections [2] O 

Briefly speaking, the symmetry method deals with the so-called canonical conservation laws 

D t pn = D x 9 n , D t p n = Dj n , n = 0, 1, 2, . . . , (1.4) 

where Dt is the evolutionary derivative and D x is the total derivative with respect to x. In 
particular, po = —F U2 /3, po = —G V2 /(3a). The recursion relations for the canonical conserved 
densities p n and p n are presented in Section [2[ All canonical conserved densities are expressed 
in terms of functions F and G. That is why equations (]l-4j) impose great restrictions on the 
forms of F and G. Equations (jl.4p are solvable in the jet space iff 

E a D tPn = 0, E a D t p n = 0, a = 1,2, n = 0,1,2,... (1.5) 

(see [28] . for example). Here 

n=0 " 

is the Euler operator. 

Conservation law with p = D x x, = Dtx is called trivial and the conserved density of the 
form p = D x x is called trivial too. This can be written in the form p € Im D x , where Im = Image. 
If p\ — P2 £ lmD x , then the densities p\ and P2 are said to be equivalent. 

There are a lot of systems in the following form 

Ut — u xxx + F(u, U x , U xx ) , Vt — a v xxx + G(u,v, u x , v x , u xx , v xx ) , 

satisfying the integrability conditions (jl.4p . Such systems containing one independent equation 
are said to be triangular. It follows from the integrability conditions that the equation for u must 
be one of the known integrable equations (KdV, mKdV etc). The second equation is usually 
linear with respect to v , v x and v xx . Triangular systems do not possess any Lax representations 
and are not integrable in this sense. Therefore triangular systems and those reducible to the 
triangular form have been omitted as trivial. 
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The system of two independent equations 

ut = u xxx + F(u,u x ,u xx ), v t = av xxx + G(v,v x ,v xx ), 

will be called disintegrated. It is obvious that the disintegrated form is a partial case of the 
triangular form. Therefore the disintegrated systems and those reducible to them have been 
omitted. 

System (jl.ip will be called reducible if it is triangular or can be reduced to triangular or 
disintegrated form. Otherwise, the system will be called irreducible. 

Our computations show that for irreducible integrable systems (jl.ip parameter a must belong 
to the following set: 

A={0, -2, -|, -| + fV5, -1-lVh). 

These values were found first in [3] and were repeated in |29j . The value of a is always defined 
at the end of computations when functions F and G have been found and only some coefficients 
are to be specified from the fifth or seventh integrability conditions (see example in Section f3 . 1 () . 
This means that it is enough to verify conditions ljl.5R for n = 0, . . . , 7 and a = 1, 2 to obtain F, G 
and a. But for absolute certainty we have verified conditions (jl.5p for n = 8, 9 and a = 1, 2 for 
each system. 

The presented set A consists of zero and two pairs (a, a -1 ). The transformation t' = at, 
u' = v, v' = u changes the parameter o ^ in (II. ip into a -1 . That is why one ought to 
consider the values {0, — |, — | + |\/5} of the parameter a. Integrable systems with a = 
were mentioned above, see also |30j . This paper is devoted to investigation of the case a = — 1/2 
only. The case a = — | + §\/5 will be presented in another paper. 

Section [2] contains recursion formulas for the canonical densities. The origin of the notion, 
some examples and a preliminary classification are considered. 

A list of forty integrable systems and an example of computations are presented in Section [3j 

Section H] contains differential substitutions that connect all systems from the list. The 
method of computations and an example are considered. It is shown that all systems from the 
list presented in Section [3] can be obtained from (jl.2p and (jl.3p by differential substitutions. 

Section [5] is devoted to zero curvature representations. The zero curvature representations 
for systems (11.21) and (jl.3p are obtained from the Drinfeld-Sokolov L-operators. A method of 
obtaining zero curvature representations for other systems is demonstrated. 



2 Canonical densities 

One of the main objects of the symmetry approach to classification of integrable equations is the 
infinite set of the canonical conserved densities. Let us demonstrate how canonical conserved 
densities can be obtained the from the asymptotic expansions for eigenfunctions of the Lax 
operators. The simplest Lax equations concerned with the KdV equation 

ut = 6uu x — u xxx 

take the following form 

*l) xx -utp- fi 2 ijj = 0, (2.1) 
ipt = ~^xxx + 6utp x + 3u x ip + 4/i 3 V>. (2.2) 

Here u is a solution of the KdV equation and [i is a parameter. The standard substitution 
tp = exp [ pdx 
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reduces equations (|2.ip and (|2.2p to the Riccati form 

p x + p 2 - u - p 2 = 0, (2.3) 
d t j pdx = -4(d x + pfp + 6up + 3u x + 4p 3 . (2.4) 



Differentiating temporal equation (|2.4p with respect to x one can rewrite it, using (|2.3p . as the 
continuity equation: 

Pt = d x [(2u-Ap 2 )p-u x }. (2.5) 
To construct an asymptotic expansion one ought to set 

oo 

p = p + Y J Pn(-2p)- n . (2.6) 

n=0 

Then equation (|2.3p results in the following well known recursion formula [2] 

/o n+ i = D x p n + ^ piPn-i, n = 1, 2, . . . , po = 0, pi = -u, (2.7) 

i=l 

and (|2.5|) results in infinite sequence of conservation laws: 

D t p n = D x (2up n - Pn+2), n>0. (2.8) 

We change here dt — > Dt and d x — > D x because u is a solution of the KdV equation. The 
obtained conservation laws are canonical. It is easy to obtain several first canonical densities: 

P2 = ~ui, p 3 = u 2 -u 2 , Pi = D x (2u 2 - u 2 ), 

It is shown in [2] that all even canonical densities are trivial. Note that if one chooses an- 
other asymptotic expansion, for example, in powers of p -1 instead of (|2.6p . then another set of 
canonical densities is obtained, which is equivalent to the previous set. 

The canonical densities that follow from (|2.7p can also be obtained by using the temporal 
equation (|2.4|) only. Indeed, setting dt J pdx = 9 one obtains from (|2.4p 



-4(d x + p) 2 p + 6up + 3u x + Ap 3 = 9. (2.9) 

Using the same expansions as above 

00 00 
P = ^ + E Pn(~W n , = 0n(-2p)- n , 

n=0 n=0 

one can obtain from (|2.9p the following recursion relation: 

n+l n 

Pn+2 = 2up n + 2 ^2 PiPn-i+1 ~ | ^ PiPjPn-i-j ~ 
i=0 i,j=0 

+ 2D X (p n+ l - ^ PiPn-i^J ~ \D\Pn ~ uSn-t + Ui5 n0 , U = -2, -1,0,..., 

where is the Kronecker delta. The obtained relation provides po = 0, p\ = —u, p 2 = 
—u\ — 9q/S, etc. As D t po = D x 9q and po = 0, then 9q = 0. The higher canonical densities p n , 
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n > 2 depend on 6 n —2' The fluxes n must be defined, now from equations (jl.4p . For example, 
9 l = u 2 - 3n 2 . 

The traditional method to obtain the canonical densities for an evolution system [25] 

u t = K(u, u x , . . . , u n ), u(t,x)£R m , m>l, u% = d k x u a . (2.10) 

consists, briefly, in the following. The main idea is to use the linearized equation 

(A--K*)V = (2.11) 

or its adjoint 

(D t + K+)cp = (2.12) 

as the temporal Lax equation. Here 

r)K a clK@ 

n,f3 du n n,f3 

The spatial Lax operator (formal symmetry) was introduced in [25] as the infinite operator 
series 

N 

R= R k°l N >^ ( 2 - 13 ) 

k=— oo 

commuting with Dt — K*. are matrix coefficients depending on u, u x , .... It was shown 
that Tr resi? (resR = R-i) is the conserved density for system (12.10j) . Canonical densities have 
been defined by the formulas 

p n = Tr res R n , n = 1, 2, . . . , 



see [26] for details. 

Operations with operator series (12.130 are not so simple, therefore we use an alternative 
method for obtaining the canonical densities. It was proposed in [32] heuristically and we 
present the following explanation (see also |33j). 

Observation. One can obtain equation (|2.9p from (|2.2p by the following substitution 



j pdx + Odt, (2.14) 



where pdx+6 dt is the smooth closed 1-form, that is, Dtp = D x 9. This implies e~ u) Dt& ) = Dt+9, 
e~ w D x e w = D x + p and so ()2.9[) follows. Another way to obtain the same equation is to prolong 
the operators Dt — > dt + 6, D x — > d x + p in (12. 2j) formally and to set tp = 1. For systems, one 
must set ip a = 1 for a fixed a only. 

We shall apply this method to system (jl.l|) now. 

The linearized system (jl.ip with prolonged operators D x — ► D a . + p, D t ^> D t + 9 takes the 
following form: 

[(A, + p) 3 + F U + F U1 (D x + p)+ F U2 (D x + pf -D t -e] ^ 
+ [F„ + (£>« + p) + F V2 (D x + p) 2 } ^ 2 = 0, 
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[G u + G U1 (D x + p) + G U2 (D x + pf] ! + a (D x + p) 3 ^ 2 

+ [G v + G V1 (D x + p) + G V2 (D x + pf - D t - 9] # 2 = 0. (2.15) 

If one sets here \&i = 1, then the first equation takes the following form 

(D x + pfp + F u + F UlP + F U2 (D X + p)p-9 

+ [F v + F Vl (D x + p) + F V2 (D x + pf] ^ 2 = 0. 

It is obvious from this equation that the following forms of the asymptotic expansions are 
acceptable: 

oo oo oo 

p = p- 1 + ^2pnv n , e = p~ 3 + ^2e n p n , m 2 = ^p n p n . 



n=0 



n=0 



n=0 



Here p is a complex parameter. Then, after some simple calculations, the following recursion 
relations are obtained (n ^ —1): 

n+l n 

Pn+2 = \&n — ^2 PiPn-i+1 — | PiPjPn-i-j ~ 3 F Ul (6 n -l + j0„) — 3 F u <5 nj 

i=0 i+j=0 



i=0 



l(F v + F Vl D x + F V2 D 2 x )ip n - \ F U2 yD x p n + 2 p n+i + ^ pip n . 

I n n 

\Fv 2 ^n+2 + 2 ^ pnp n -i + i + ^ PiPjVn- 



-1-3 



i=0 



i+j=0 



\F V2 2 (/9 n+ i + ^ pjDj; </? n _j + D x Pifn-i 



i=0 



i=0 



As 



i=0 



Pn+1 + 3 Ac Pn + | PtPn-i 
i=0 

n 

(1 - a)ip n+3 = G u 5 n>0 + G U2 (D x p n + 2 p n+i + ^ Pip n -%) + G Ul (£ n -i + Pn) 

i=0 

n / n \ 

_ X/ ^n-i + G„ <£ n + G Wl Ac^n + Vn+l + ^ Pi^n-i J "Aft 
i=0 V i=0 / 

(n n \ 

2 AcVn+l + X] ^* ^ ^ n ~« + ^ J 
i=0 i=0 / 

/ n+l n 

Pitfn-i+l + PiPjPn- 



-i-J 



i=0 



i+j=0 



n+l 



n+2 



+ a A^<Pn + 3 a A^<Pn+i + 6 a Pi Ac <p n -i+i + 3 a ^ pnp n - i+2 



i=0 



i=0 



+ 3 a D x cp n+2 + 3 a Ac Vn-t-j + 3 a <Pn-i+l Ac Pi 

i+j=0 i=0 
n n+l 

+ 1° Vn-i-j D x (piPj) + 3a ^ pipjip n -i- j+1 
i+j=0 i+j=0 
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PiPjPk fn—i—j—k- 

i=0 i=0 i+j+k=0 

Here 5i k is the Kronecker delta, F Ul = dF/du\ and so on. From the recursion relations it is 
obvious why the value a = 1 is singular. Some of initial elements of the sequence {p n , (p n } read 

Po = —kF U2 , (^0 = 0, <f i = G U2 , 

° 1 — a 

others are introduced via the 5-symbols. 

If one sets in (|2.15p ^2 = 1 and a^O, then one more pair of recursion relations for {p n , (f>n\ 
is obtained. These recursion relations give us any desired number of canonical densities. As an 
example, we present here some more canonical densities: 

1 1 2 1 11 

Po — ~~ g F U2 , pi — — F U2 — — F Ul + 7^ F V2 G U2 + — D x F U2 , 

Po = — 7T~ G V2 , Pi = Tr~n G 2 — —— G V1 — - — - F V2 G U2 + — D X G V2 , (2.16) 

6 a ya z da Sab Sa 

where b = a — 1. The tilde denotes another sequence of canonical densities. Further canonical 
densities are too cumbersome, therefore we do not present them here. 

To simplify investigation of the integrability conditions, an additional requirement is always 
imposed. This is the existence of a formal conservation law \25\ [26] . A formal conservation law 
is an operator series in powers of D~ x . An equation for the formal conservation law can be 
written in the following operator form 

(D t - K*)N = M{D t + K* + ). (2.17) 



The form of this equation coincides with the form of the equation for the Noether operator 
That is a formal conservation law may be called a formal Noether operator. 

If (Dt — K*,L) is the Lax pair for an equation, then (Dt + A'* + ,L + ) is obviously the Lax 
pair for the same equation. Hence, canonical densities obtained from (|2,lip must be equivalent 
to canonical densities obtained from (|2.12p . 

It was shown in [21] that the first sequence of the canonical densities p n for system (jl.ip 
obtained from (12. lip is equivalent to the first sequence of the canonical densities r n obtained 
from (|2.12p and the second sequence of the canonical densities p n is equivalent to the second 
sequence of the canonical densities f n . Hence, p n — r n G IulL^ and p n — f n G IraD x , or 

E a (p n -r n )=0, E a (p n -f n )=0, a = 1,2, n = 0,1,2,.... (2.18) 

Equations (|1.4p (or (II. 5p ) and (I2.18P are said to be the necessary conditions of integrability. We 
shall refer to it simply as the integrability conditions for brevity. 
Our computations have shown that 

to = -Po, 7b = -p , ti = Pi, ri = pi. (2.19) 

Other "adjoint" canonical densities Tj and essentially differ from the "main" canonical den- 
sities pi and p~k- All canonical densities can be obtained using the Maple routines cd and acd 
from the package JET (see |36j). These routines generate the "main" and the "adjoint" canon- 
ical densities, correspondingly, for almost any evolutionary system (an exclusion is the case of 
multiple roots of the main matrix of the system under consideration) . 

Thus, according to (12.161) and (|2.19p we have F U2 G ImZ^ and G V2 G ImZ^ (a / 0). This 
implies the following lemma. 
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Lemma 1. System (jl.ip with a(a — 1) ^ satisfying the zeroth integrability conditions ()2. 18|) 
reads 

3 3 

Ut = u 3 - — u 2 D x f + — f Ul ul + F 1 (u,v,u 1 ,vi,v 2 ), 

v t = av 3 - -^-v 2 D x g + g Vl v% + G 1 (u,v,u 1 ,vi,u 2 )- (2.20) 
25 4g 

where ord(/,g) ^ 1. 

Indeed, one may set F U2 = —3/2D x In / and G V2 = —3/2aD x Ing, where ord(/, g) ^ 1 because 
ord(F, G) < 2. Then equations (I2T2H follow. 

From higher integrability conditions one more lemma follows. 

Lemma 2. Suppose system (|2.2U|) is irreducible and satisfies the following eight integrability 
conditions p 2 — r 2 S \m.D x , p 2 — f 2 E \m.D x and D t p n G ImD x , D t p n £ ImD I; where n = 1,3,5. 
T/ien the system must have the following form 

3 3 
ut = u 3 - — u 2 D x f + — f ui ul +fiv% + h V2 + h, 

v t = av 3 -^v 2 D x g + ^g Vl V2+giul+g 2 u 2 + g 3 , o^O, (2.21) 
25 % 

where ord(f,g,fi,gj) ^ 1. 

A scheme of the proof has been presented in [21] . 



3 List of integrable systems 

As it is shown in Section [2] the problem of the classification of integrable systems (|1.1|) is 
reduced to investigation of system (|2.2ip . That is why it is necessary to start by investigating 
its symmetry properties. 

Lemma 3. System (|2.2ip are invariant under any point transformation of the form 

(a) t' = a 3 t + P, x 1 = ax + jt + 5, it' = u, v' = v, 

(6) u' = hi(u), v' = h 2 (v), 

and under the following permutation transformation 

(c) t = at, u = v, v = u, 

where a, /3, 7 and 5 are constants, hi are arbitrary smooth functions. 

The classification of systems of type (|2.2ip has been performed by modulo of the presented 
transformations. 

Moreover, some systems (|2.2ip admit invertible contact transformations. An effective tool 
for searching such contact transformations is investigation of the canonical conserved densities. 
For example, system (|3.24p from the next section has the first canonical conserved density of 
the following form: 

Pl = { Vl -lue v ) 2 + 2c\e- 2v . 

It is obvious that the best variables for that system are 

U = e~ v and V = «i - \ue v . 
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This is an invertible contact transformation. In terms of U and V the system takes the following 
simple form: 

U t = D x (U 2 + \UVy - \UV 2 + \c\U 3 ) , 
V t = \D X (V 3 - 2V 2 ) - \c\D x (WXJ x + U 2 V). 

If ci ^ this system can be reduced to f|3. lQj) by scaling, otherwise the system is triangular: the 
equation for V will be independent single mKdV. Moreover, the equation for U becomes linear. 
That is why a / in (pT24]) . 

Canonical densities for the triangular systems contain only one highest order term in the 
second power as in the considered example p = V 2 or p = V 2 + • • • , or p = V xx + • • • etc. Trian- 
gular systems and those reducible to the triangular form have been omitted in the classification 
process as trivial. 

To classify integrable systems with a(a — 1) 7^ one must solve a huge number of 

large overdetermined partial differential systems for eight unknown functions of four variables. 
This work has required powerful computers and has taken about six years. All the calculations 
have been performed in the interactive mode of operation because automatic solving of large 
systems of partial differential equations is still impossible. The package pdsolve from the 
excellent system Maple makes errors solving some single partial differential equations. The 
package diffalg cannot operate with large systems because its algorithms are too cumbersome. 
Thus, one has to solve complicated problems in the interactive mode. Hence, to obtain a true 
solution one must enter true data! Under such circumstances errors are probable. The longer 
the computations the more probable are errors. This is the reason why we cannot state with 
confidence that all computations have been precise all these six years. That is why the statement 
on completeness of the obtained set of integrable systems is formulated as a hypothesis. 

In this and in the following sections c, c,, k, ki are arbitrary constants. 

Hypothesis. Suppose system (12.21|) with a = —1/2 is irreducible. If the system has infinitely 
many canonical conservation laws, then it can be reduced by an appropriate point transformation 



to one of the following systems: 

u t = u 3 + vui, v t = — 5^3 + uui — vv\; (3.1 

ut = u 3 + v± ui, v t = -\v 3 + \(u 2 - uf); (3.2 

Ut = U3 + vui, v t = — — vvi + ui; (3.3 

ut = «3 + !;Mi+t;i u, v t = - vv\ + u; (3.4 

u t = u 3 + v 1 ui, v t = -\v 3 - \v\ + u; (3.5 

Ut = 113 + uui + vi, v t = -\v 3 + |«iti2 - uvi; (3.6 

ut = U3 + v 2 + kui, v t = — \vz + |un 2 + \u\ + \v? + k [u 2 - v\) ; (3.7 

ut = u-s + 2%i>2 + ft>i + 3V - k [v +ui), v t = -2V3 + U2 + kvi; (3.8 

u t = u 3 - ^U!V 2 - juxvf + \u\, v t = -\v 3 + |uiu 2 - \u\vi + \v\ \ (3.9 

u t = {u 2 - I uvi - \uv 2 + ju 3 ) x , v t = {-\v 2 + §mii - \u 2 v + \v 3 ) x ; (3.10 

ut = u 3 - \v 2 - \u\v x - \u\, v t = -\ v 3 + I (vi - u 2 + \u\) 2 - I vj; (3.11 
u t = {u2-\vi-\uv -\u 3 ) x , v t = (-\v 2 + l{v -ui + \u 2 ) 2 -\v 2 ^; (3.12 
u t = u 3 - 3gv 2 - 3ui(m + vt) - | v\ - <ov x g 2 - c x g 3 - 3g 4 , 

vt = ~\v 3 - jCiu 2 + 3uj - tjV 2 - 6ui5 2 + c 1 g 3 + 3g A , g = u + v; (3.13 

Ut = u 3 — 3u\v\ + (u — 3v 2 )m, vt = — \v 3 + \ii2 — u\v — (u — 3v 2 )v\; (3.14 

ut = u 3 — 3u\v 2 + uu\ — 3uivf, vt = — \v 3 + \u\ — u v\ + vf; (3.15 
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ut = 1*3 — | u\V2 — | mm + n ui — c me - ' 



ut = «3 + 3ii2Ui + |uif2 + I uiv\ — uu\e 2v — e c 
v t = 

ut = u 3 + 3u 2 vi + |uif2 + fuiu 2 - uu\e 2v - \u\e~ 2v 



ut = u 3 + (k + \J u 1 + Vi j Ui, 

vt = 4 « 3 + | (2 "y, + " 2)2 " 3uu 2 - A;(2m 2 + Wl ) - |(n 2 + m) 3/2 ; (3.16) 
u -\- V\ 

n(2vV\+U2) 2 q 9 o q , . 9 v 

«t = «3 ~ Z " 5 — + 3w 2 + 4«i + 4^ - M 2 ^ + Ml , 

4 IT + Ml 

u t = 4^3 + 5 U 2 + kvr, (3-17) 
Ut = ""3 + Ul ^+vf ~ t UlUl + ciUiy/u+Vi, 

Vt = 4«3 - |U 2 + | (m+^) 2 + 2^2 _ 4 M 2 _ 2ni ^ + ^ _ | Cl ( M + ^3/2. (3 . 18) 

ut = us + iti \A* + ^i — kui, 

v t = -i«3 - §«2 + | ( ^ + ^ )2 - |(u + ui) 3 / 2 + 2A;n + fet, i; (3.19) 

ut = ""3 + uvi + (u 2 + v)ui, v t = -\v% + 3iiiii 2 - (u 2 + v)v±; (3.20) 
ut = u 3 + 3(u + k)v 2 + 3iii(t>i + u 2 ), 

vt = - j«3 - §uii2 - §(ui + u 2 ) 2 - \u\ + /cu 3 + § u 4 ; (3.21) 
ii t = u 3 — §t> 2 - | uivi - \u\ - 3ui(cie u + 2c 2 e 2lt ), 
u t u = -i«a + | (\u\ -u 2 + v 1 + ae u + 2c 2 e 2 ' 11 ) 2 

-\v\-\c x u2e u + \c\e 2u + 2 Cl c 2 e 3u , a + or c 2 + 0; (3.22) 



u* = 4«3 + i u i + Ul ~ U Vl + c2yi e ~ 2v > ( 3 - 23 ) 

^t = us — §mu 2 — | ttif 2 + ui e"(ui + 2uv\) — jr u 2 uie 2i; — |c 2 ui e~ 2v , 

v t = -\vz + v ^ + u 2 e v + \ue 2v {2u 1 +uv 1 ) + lc\v 1 e- 2v , a ^ 0; (3.24) 

,3 ,9 2 2?) -3« 



-^ 3 + 3 u 3 + (ni+-ut;i)e 2,; ; (3.25) 



ut = 4u 3 + \v\ + ( Ul + m^e 2 " + iuie" 2 "; (3.26) 
ut = u 3 + \u\V2 + 3n 2 ui + \u\v\ - c\ u\e~ 2v - \ u\e 2v {u 2 + c 2 ), 

v t = 4u 3 + 1^1 + c\v ie - 2v + \ e 2v {2uui + tt 2 ui + c 2 ui); (3.27) 
ut = u 3 + §mu 2 + 3n 2 ui + |uiu 2 - i e 2v (ui(6u 2 + ci) + 4nvi(2n 2 + ci)) 
+ (v 2 (2u 2 + ci) + (ui + 2nvi) 2 + 2ci?; 2 ) , 

u f = 4«3 + i u ? + 3 e<2v ( 4u< "i + ( 6n2 + c i) v i) + ^ («2 + 2n i y i) ; ( 3 - 28 ) 

ut = u 3 + | uiu 2 + 3u 2 ui + | uiu 2 + 3^^^'2(cl^te , ' + c 2 ) + ci(c 2 - l)n 4 e 3v 

- | u 2 e 2v {ui{\ + 5c 2 ) + 8c 2 wi + 2c 2 u(l - 3c 2 )) - 3c|(ui + 2u«i) 
+ ^c\e v {ui + 2itui — 2c2u) 2 + |c2Ui(2ni + 3«ui), 

ut = 4u 3 + 2 U ? + \c\e"{u 2 + 2mui) + c\{\ - cj)u 3 e 3v + 6cic 2 ue ,; (7;i - c 2 ) 

+ f ue 2t, (2ni(l + c\) + uu^l + 5c 2 ) + 2c 2 u(l - 3c 2 )) + |c 2 ui(2c 2 - ui); (3.29) 
ut = u 3 + I uiu 2 + 3-u 2 ui + | Miu 2 + 3e v (u 2 + c)(u 2 + 2v\) + | e^tti(ui + Auv\) 

- \e 2v ((3n 2 + c)ui + 4(n 2 + c)u«i) , 

ut = -|u 3 + \ v\ + § e w (^ 2 + 2mm) + |e 2,) (2«m + (3u 2 + c)ui) ; (3.30) 
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u t = u 3 - - \u\v\ + \u\ -ae 2v u x - c 2 (tti + 2vi) e 2{u+v) 

+ c 3 { Ul -2vi)e 2{v - u \ 
vt = ~\v 3 + %U X U2 - |u?t;i + \v\ + d e~ 2v Vl + (c 2 e 2 ^ - c 3 e 2 ^) Vl ; (3.31) 
u t = u 3 - ^u\v 2 - juivf + \u\ + (c 2 e" + c 3 e~ u - 3c? e~ 2v ) u±, 
v t = -\v 3 + jjmua - + j«i + (c 2 - c 3 e~") m 

+ (3c? e~ 2v -c 2 e u - c 3 e" u ) «i ; (3.32) 
u t = u 3 - %U!V 2 - + \u\ + 3 fc(u? - 2v 2 ) e~ v - 3(c 2 - 3k 2 ) Ul e~ 2v 

+ 8k(k 2 -c 2 )e- 3v , 

v t = ~\v 3 + \u x u2 - \u\vi + \v\ + 3k(u 2 - 2u lVl ) e~ v + 3(c 2 - 3k 2 ) Vl e~ 2v ; (3.33) 
u t = u 3 - \u x v 2 - \u x v\ + \u\ - 3cj Ul e 2( - u+ ^ + 3du 1 {u 1 + 2v 1 )e u+v 

+ (c 2 e-"-3c^ e- 2v ) Ul , 
vt = ~\v 3 + lu x u 2 - lu 2 lVl + \v\ + 3c 2 (2m + «i) e 2 ^ + 3ci(u 2 + uj)e u+v 

-c 2 {u l +v l )e~ u + 34v ie - 2v - (3.34) 
u t = u 3 - luiv 2 - \uivl + \u\ + 3c 2 ui{ Ul + 2vi) e u+v - 4cic 2 e^ u+v) 

+ 3[{c l -c 2 2 )u l + 2dv 1 }e 2 ^ +v \ 
vt = ~\v 3 + \uiu 2 - \u\vi + \v\ + 3c 2 (u 2 + u?) e u+,; + 4cic 2 e 3{u+v) 

+ 3[2c£«i - (ci - cl)^i] e 2(n+,;) ; (3.35) 
u t = u 3 - \u\v 2 - \u\v\ + \u\ + I c\u\ e~ 2v + d{2v 2 - u\) e~ v 

- 2dc 2 (u 1 + 2vi) e u + 3c 2 Ul ( Ul + 2«i) e n+,; - 3c^ m e 2 ( M +^, 

i>i = + |«i«2 - + 3^1 - § c?wi e -21 ' + d{2u\vi - u 2 ) e~ v 

- 2dc 2 (ui - v x ) e u + 3c 2 (u 2 + u\) e u+v + 3c\ (2ui + Vl ) e 2{u+v) ; (3.36) 
n t = u 3 - \ Ux v 2 - \ u lV j + \ u\ - 3 Ul [c? e 2 ^ + c 2 e 2 ^ + 2 Cl c 2 ' 

- 3c§ ui e~ 2t ' + 3ci«i(«i + 2ui) e u+1 ' - 3^ui(«i - 2«i) e^ - ", 
Ut = -5U3 + |«ii*2 - |^ 2 ui + i^f + 3c 2 (2ui + vi) e 2(u+,;) + 6cic 2 ui e 

+ 3c| (ui - 2ui) e 2iv ~ u) +3d(u 2 + n 2 ) e u+v + 3c 2 (n 2 - u 2 ) e v ~ u + 3c§ «i e~ 2v ; (3.37) 
«t = us - |«i«2 - + |«? - 6c? e 3 ^) - I c 2 (5m - 8«i) e 2 ^) 

+ § dui(ui + 2«i) + 2c? c 2 e 2u+,; + § cic| e""" - ± cic 2 (7ni + 12«i) e u 

+ c 2 (2v 2 - u?) e -« + 11 c 2 m e" 211 - | c 3 e" 3 ", 
^ = -1^ + | Ul u 2 - fufo + l^f + 6c? e 3 ^) + I c?(18ni + 5«i) e 2 ( u+,; ) 

+ § d(u 2 + uf) e u+v - 4c? c 2 e 2 ^ + \cy<% e u ~ v - \ dc 2 (u x - Vl )e u 

+ c 2 (2niui - u 2 ) e~ v -%c\v x e~ 2v ; (3.38) 
u t = u 3 - lu lV2 - \u x vl + \u\ + c 3 e- v {u\ - 2v 2 ) + \c\u x e- 2v 

+ ci (3 Ul e u+v + 2c 3 e u ) (ui + 2ui) - c 2 (3nie ,J - ti + 2c 3 e- M )(-ui - 2^) 

-3m(cie n + c 2 e- n ) 2 e 2 ", 
ut = -5^3 + |uiii 2 - |«f«i + lyf + c 3 e~ v (u 2 - 2mvi) - 2 l clv 1 e~ 2v 

+ 3c?e 2( ' il+1,) (2ui + vi) - 3c 2 2 e 2 ^- u \2 Ul - v x ) + 3cie u+ > 2 + u \) 

- 3c 2 e v ~ u (u 2 - uf) + 2c 2 c 3 e""(m + «i) + 2cic 3 e n (tii - v x ) + 6dc 2 vie 2v ; 



(3.39) 
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3 (2ff 3 - u 2 + 2gvi \ £ 2 q 3 2 

iH = u 3 - f 5 h 6g(u 2 - v 2 ) - 6^ - 9ui?;i - 

- 6{5g 2 + 4cg + c 2 )ui - 6g 2 vi + 2cg 2 (8g + 3c) + 9g 4 , 

l . 3 ( 2 5 3 - w 2 + 2^i) 2 3 2 2 2 ^ 

«t = + § 5 3(35 + c)u 2 - f^i + 6{9g + 8cg + 2c )u\ 

u\ — g z z 

+ 6(6g 2 + 4cg + c 2 )v 1 -2cg 2 (8g + 6c) -9g 4 , g = u + v; (3.40) 

Remark 1. Systems (]3.ip . fl3.3f) . (j3.6f) and (|3.20j) were proposed in [3], where system (|3.20p is 
given with a misprint. System (|3.10p was presented in [29j . 

Remark 2. Ten pairs of integrability conditions (for po~P9 and po~pg) have been verified for 
each system (|3.ip - (|3.40p . and nontrivial higher conserved densities with orders 2, 3, 4 and 5 
have been found. 

Remark 3. It is shown in [21] that systems (|3.10p and (|3.12p are unique divergent systems of 
the form (|2.2ip that satisfy the integrability conditions. 

Remark 4. System (I3.22p is a modification of (|3.1ip . systems (|3.3ip - (l3,39p are modifications 

of dS^D . 

Remark 5. Canonical densities for system (13.251) depend on the nonlocal variable w = D~ 1 e~ v . 

Remark 6. Many of the systems possess discrete symmetries. They are: 
u -> -u for §n$), (I3TTU1) . (pngp , (IQl and (ED ; 

it — ► —u, v — > v -\- iri for (|3.24p ; 
u —> iu, v — > f — |tt, ci — > — ci for ()3.28|) : 

{u — ► — tt, f — > v + iri} U {v — ► v + wi, c\ — > — ci} U {u — > —u, c\ — > — ci} for (j3.29H : 

{u — > —u, v — > v + vri} U {u — > iu, v — > v — |7r} for (|3,30p ; 

u — > -u, c 2 — > c 3 , c 3 — ► ci for (|3.32j) ; 

-u — > —it, A; — > — /c for (|3. 33|) ; 

it — > —it, ci — > C2, C2 — > ci for (|3.37|) ; 

it — > —u, C3 — > -C3, C2 — > ci, c\ — > C2 for (|3.39p . 

Also, systems f)3. 1 [) . (j3.2|) . (|3.9p . f|3. 10|) and (|3.27p preserve the real shape under the transfor- 
mation u — > iu. System f|3.33j) keeps the real shape under the transformation u — > iu, k —> ik. 

3.1 Example of computations 

Let us consider the simplest case of system (jl.ip : 

Ut = u 3 + fi(u,v)ui + / 2 (u, v)vi, v t = av 3 + g 1 (u,v)u 1 + g 2 (u,v)v 1 , (3.41) 

where a(a — 1) 7^ 0. Formulas (|2.16p are reduced now to the following 

11 , 

Po = 0, po = 0, px = --f 1 , pi = - —92- (3.42) 

3 6a 

The further canonical densities read 

P2 = + /2, u «l) + r Ar/l, p 2 = - — (gi, v ui + 92, v vi) + ^-D x g 2 , 

3 3 6a 6a 

T2 = -(/l/a^l + f2,uV\), f 2 = — {g\,vU\ + 52,v«l), (3.43) 

where indices after commas denote derivatives. 
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The first integrability condition (|1.5[1 for p\ can be split with respect to U3, v 3 , u 2 and v 2 . 
This provides the following equations 

fi,uv = 0, fi )U uu = 0, fi >VW) = 0, or fi(u, v) = ciu 2 + c 2 u + c 3 v 2 + c 4 t> + c 5 . 

Analogously, condition (|1.5p for pi implies 32 (u, v) = b\u 2 + b 2 u + 63V 2 + 64U + 65. It is obvious 
from (|3.43p that the second integrability conditions (|2.18p are r 2 £ ImD x and f 2 £ lmD x . These 
conditions provid e f 2 , u u = 9i,w = or f 2 (u, v) = uf 3 (v) + / 4 (v ), 51 (u, u) = 1/03 (tt) + g${u). 
Thus, system (I3.4ip takes the following form: 

u t = u 3 + (ciu 2 + c 2 u + c 3 v 2 + c 4 v + c 5 )«i + (tt/ 3 (v) + fi{v))vi, 

v t = av 3 + (vg 3 (u) + g±(u))ux + (b\u 2 + 6 2 u + 03t> 2 + hv + 05)^1 • (3.44) 

Now one can obtain 6 2 = D~ 1 Dtp 2 and 6 2 = D~ 1 Dtp~ 2 in an explicit form. The expres- 
sions D t p\ and D t p\ are not the total derivatives yet: 

D t p\ = D x hi(ui,Vj) + Rt(u,v,ui,vi) = D x 0i, 
D t pi = D x hi(ui,Vj) + Rt(u,v,ui,vi) = D X Q X - 

Therefore, we have set 9\ = h\(ui,Vj) + qi(u,v), 9\ = h\(ui,Vj) + qi(u,v), where q± and q± are 
unknown functions and D x q\ = Ri, D x q\ = R\. This trick allows us to evaluate r 4 , /3 4 , T4 
and verify the fourth integrability conditions f|2. 18j) . These conditions imply /g = g 3 = 0, hence 

f 3 = aiv + a 2 , g 3 = a 3 u + a 4: . 

To simplify the further analysis one must list all irreducible cases of f\ (or g 2 ). Let us take 
fi = c\u 2 + c 2 u + c 3 v 2 + c^v + C5 for definiteness. 

Lemma 4. Using complex dilatations of u and v, translations u — > u + Ai, v — > i> + A2 and i/ie 
Galilei transformation ut — ► ih + Q^i, ^ — ► nj + at^ one can reduce f\ to one of the following 
forms: 

l)u 2 + v 2 ; 2)u 2 + av; 3) v 2 + au; 4) u + v; 5) tt; 6) n; 7) /1 = 0, 

where a is any constant. Moreover, in the cases 4-7 £/te function g 2 must be linear (61 = 63 = 0) 
because otherwise the permutation u <-> v gives one of the cases 1-3. 

In cases 1 and 3 contradictions follow from the integrability conditions (jl.5p with n = 1,3 
and (|2.18p with n = 2,4. In case 2 the integrability conditions (|1.5p with n = 1,3,5 and (|2.18p 
with n = 2,4 are satisfied iff system ()3.44p is reduced to a pair of independent equations. Thus, 
a nontrivial integrable system (|3.41|) must belong to the following class: 

Ut = u 3 + (c 2 u + Civ)ui + (u{a\v + a 2 ) + f±{v))v\, 

v t = av 3 + (v(a 3 u + 04) + g±{u))u\ + (b 2 u + & 4 i> + b 5 )v±, (3.45) 

and only the following cases are possible: 

4) c 2 = c 4 = 1; 5) c 2 = 1, c 4 = 0; 6) c 2 = 0, c 4 = 1; 7) c 2 = c 4 = 0. 

In case 4 the integrability conditions (jl.5p with n = 1, . . . , 5 provide the functions 54 = k\u + k 2 , 
/ 4 = &3t> + fc 4 , the coefficients 03 = 0, a 4 = — 1 + a 2 + b 2 , 04 = 6 2 = (a + l)o 2 — 2a — 1 and the 
following equations: 

(a + l)(2a 2 + aa 2 - 2a - 1) = 0, a 2 (a 2 - 2) - a(4a| - 5 - 7a 2 ) + 16a 2 - \\a\ -3 = 0, 
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a 5 (a 2 - 2) - a 4 (4a| - 7a 2 - 6) - 2a 3 (16a 2 , - 5a 2 - 37) 

- a 2 (177a?, - 224a 2 - 41) - a(236a| - 353a 2 + 102) - a 2 (37a 2 - 53) - 17 = 0. 

Using the package Groebner in Maple, one can obtain a 2 = (1 — a)/3, a 2 + 7a + 1 = or 
a = (3c — 7)/2, c 2 = 5. Then, the remaining coefficients are also determined and we obtain 

ut = u 3 + {u + v)u\ + \ ((3 — c)u + (5c — ll)v)v\, c 2 = 5, 

v t = ±v 3 (3c - 7) + ((c + 2)u -v)ui + |(c - 3)(u + v)vi. (3.46) 
The following substitution 

V = \{c + l){v -u), U =±{c + 3)u + \v 
reduces system (|3.46p to the third Drinfeld-Sokolov system 



U t = -8U 3 + 3U 3 + Q(V - 8U)Ux + 12UV U 
V t = 12U 3 - 2V 3 + 48VUi + 12(2*7 - V)Vi 

that has been presented first in [3]. Scaling 

t^-\t, U=^±U, V=^-\V 

gives more symmetric form of system (|3.46aD 

U t = AU 3 + 3V 3 + (4U + V)U! + 2UV h 
V t = 3U 3 + V 3 - AVUi - 2(V + U)Vi 



(3.46a) 



(3.46b) 



that was found in [13] . 

In case 5 the equations a± = a 2 = 0, g'^f^ = 0, g'" = follow from the integrability condi- 
tions (jl.5p with n = 1, . . . , 5. This implies f^^O because otherwise the first equation of (|3.45[) 
will be independent. Hence, there are two branches (1) / 4 = 1 and (2) f' A ^ 0, g'± = 0. Along 
the first branch, if one use additionally the integrability conditions (jl.5p with n = 7 and solves 
a large polynomial system for constants, one can obtain the following system: 

ut = U3 + uux +vi, v t = -2v 3 - UVi, 

that can be transformed to (II. 3p by a scaling. 

Along the second branch the integrability conditions (jl.5jl with n = 1,...,5 provide the 
following system 

ut = u 3 + uui - vvi, v t = -2v 3 - UVl, 

that is equivalent to (|1.2p . 

Case 6 is symmetric to case 5: one can obtain f%' = 0, g'l = 0, 03 = 0,4 = 6 2 = from the 
integrability conditions (jl.5j) with n = 1, . . . , 5. Hence 54 7^ and we have two branches 54 = 1 
° r 94 — ^- Using the additional integrability conditions (|1.5p with n = 6, 7 one can obtain 
equations (11. 2p and (II. 3p . 

There are many branches in case 7 but all of them provide linear or triangular systems only. 

As one can see, classification of integrable systems of the form (|3.4ip is a sufficiently laborious 
task. System (I2.2ip contains eight unknown functions depending on four variables, therefore 
classification of these systems is much more difficult. 
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4 Differential substitutions 

A differential substitution is a pair of equations 

u = f(U,V,U x ,V x ,...,U n ,V n ), v = g(U,V,U x ,V x ,...,U n ,V n ), (4.1) 
where / and g are some smooth functions. 

Definition 1. If for any solution (U, V) of a system (£) formulas (|4.ip provide a solution (u, v) 
of system (jl.ip . then one says that system (jl.ip admits substitution (|4.ip . 

In all cases that we know, the new systems (£) belong to the same class (jl.ip 

U t = U xxx + P(U, V, U x , V x , U xx , V xx ), V t = a V xxx + Q(U, V, U X ,V X , U xx , V xx ), (S) 

with some smooth functions P and Q. There exist some group-theoretical explanation of this 
fact for KdV type equations [35J. Our attempts to introduce another parameter a' 7^ a in ([21) 
had no success. 

Substituting (|4,ip into (jl.ip one obtains the following equations 
(D 3 J + F(f, g, D x f, D x g, D 2 J, D 2 x g) - d t /)| ^ = 0, 

{aD x g + G(f,g,D x f,D x g,D 2 x f,D 2 x g)-d t g)\ =0. (4.2) 



It is obvious that transition to the manifold (S) in (|4.2p is equivalent to a replacement of <9j by 
the evolutionary differentiation Dt performed in accordance with (j^g|) : 

D t f = Dlf + F(f, g, D x f, D x g, D 2 J, D 2 x g), 

D t g = aDlg + G(f, g, D x f, D x g, D 2 x f, D 2 x g), (4.3) 

where 

D *f = E w Dl(Us + p) + E w Dl{aVi + Q) - 
1=1 1 1=1 1 

Another way to obtain (j4.3p is to differentiate equations (j4.ip with respect to t in accordance 
with (jl.ip and ((SI) and exclude u and v by using (j4.ip . This algorithm and many others are 
coded in Maple, see for example [36] . 

To find the admissible functions /, g, P, Q from (|4.3p one can use the following easily provable 
formula: 







D m f _^(™\ D m-s 9f 1L= for i >0 



dU k 

h s=0 

and the analogous formula for d/dV^. Differentiating (j4.3p with respect to U n+ z and V n+ 3, one 
obtains 

df 0, #=0. (4.4) 



Other corollaries of (j4.3p are too cumbersome to consider them in the general form. 

Let us consider, as an example, the first order differential substitutions for system (jl.2 



ut = u 3 + v ui , v t = -lv 3 + uui - vv\. (4.5) 
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According to (|4.4p one has f = f(U, V, Ux), g = g(U, V, Vx), hence equations (|4.3p now read 
Dlf - fu(U 3 + P) - fv(Q ~ V 3 /2) - f Vl (Ui + D X P) +gD x f = 0, 

gu(U 3 + P)+g v (Q- V 3 /2) + g Vl (D x Q - V 4 /2) + \ D 3 x g - f D x f + g D x g = 0. (4.6) 
Differentiating (|4.6p with respect to U 3 and V 3 one can obtain four equations: 

df dP 

dUi dU 2 
dg dQ 

dVi dU 2 

Let us consider some corollaries of these equations. 

1. If df /dU\ = and dg/dVx = 0, then u = f(U), v = g(V) is a trivial point transformation. 

2. If df/dUi = 0, then u = f(U) and one can set f(U) = U by modulo of the point 
transformation. In this case P = gU\ from the first of equations (|4.6p . 

3. If dg/dVi = 0, then v = g(V) and one can set g(V) = V by modulo of the point transfor- 
mation. In this case Q = fD x f — VV\ from the second of equations (|4.6[) . 

4. If (d f / dU\)(dg / dV\) ^ 0, then one can find P and Q as polynomials of U 2 and V 2 from 
equations (|4.7p . 

Investigation of cases 2—4 provides seven nontrivial solutions of equations (|4.6|) (see below 

(1311) -> (1321), iU) (I33ZD)- 

Note that integrable system (I3.6P admits strange differential substitutions that generate non- 
integrable systems. For example, system (|3.6p admits the following differential substitution: 

u = \V 2 - - \U x e v , 

v = I {-V A + VtVs + V?V 2 - \V? - Ufe 2V + e v (U 3 + 2U 2 V X + W X V 2 )) , (4.8) 

so that the functions U and V satisfy the following system: 

U t = U 3 + \U 2 V X + piV? - U ie - y f"(U) + /(Z7), 
Vt = -|F 3 + i^ 3 + \e v (U 2 + U X V X ) - f'(U) 

with arbitrary function /. This system does not satisfy the integrability conditions (jl.4p . To 
comprehend this unusual phenomenon we evaluate V 2 , V 3 and V4 from the first equation (14. 8p 

y 2 = lu + U x e v + \V?, V 3 = lux + D x (U ie v ) + VM, 
Vi = \u 2 + D 2 x (U ie v ) + D x (ViV 2 ), 

and substitute them into the second one. The result is 

u = — u — %u 2 . 

It is easily verified that the obtained constraint is a reduction of system (I3.6P into the single 
KdV equation ut = —l/2u 3 — uu\. This means that using a substitution like (|4.8p we are trying 
to construct an integrable system from the single KdV equation. There are other such examples 
for system (I3.6p . Note that the reduction obtained above follows from the reduction u = const 
for system (see ([33]) and -> Ijg^jl ). 



3D 



9/ 



dUx 

38g_ 
2dU' 



df dP 

dUx dV 2 
dg dQ 

WxW 2 ~ 



3 5/ 

2 ay 

2 x dV{ 



(4.7) 



Integrable Evolutionary Systems and Their Differential Substitutions 



17 



To organize the presented list of systems we have computed admissible differential substitu- 
tions for each system and present the results in this section. The formula 



u 



f(u,v,u x ,v x , ...), 



g(u,v,u x ,v x , ...) 



(A) - (B) 



will denote that if u' and v' are substituted into system (A), then system (B) follows for u and v. 
We say in this case that system (B) is obtained from system (A) by the differential substitution. 



Substitution (A) — > (B) establishes an interrelation between the sets of solutions of sys- 
tems (A) and (B): (u,v) \— > (u',v') is a single valued map. And conversely, if for some solu- 



tion (V, v') of system (A) one solves the system of two ordinary differential equations (A) — > (B) 
for u and v, then one or more solutions of system (B) are obtained. Of course, explicit solutions 
can be obtained very rarely when the substitution is linear or invertible (see below). 
Let us consider the following simple example: 



u = ui, v = V\. 
This substitution is possible for any divergent system 

ut = {u 2 + F(u,v,ui,vi)) x , v t = {~\v 2 + G(u,v,ui,v 1 )) x . 



It produces the system u t = 113 + F(ux, v±, U2, V2), vt = —V2/2 + G(u±, v±, u 2 , v 2 ) without uq 
and vq. The inverse transformation is quasi-local u = D~ l u\ v = D~ l v' . This is a well known 
fact, that is why the substitutions (u,v) — > (ui,vi) are not written for the divergent systems 
below. In some cases analogous substitutions are not so obvious and we present them likewise 
dO -> ([321), for example. 



each other. 

The proof can be obtained by a direct verification. 
List of the substitutions: 



u, v = «i; 
^= (u 2 - UlVl), 
(ux - uv), 



W2-|(«l+«l); 

2\. 



3/„,2 



- |(u + V 



V2(u 



3 



3V2((u + v 
\/2u, v 



v 
1 



2u 2 + 2ui) , 
2«i + 2v) , 

Ul ) + ±V2cl 
-3vi + u — 3v 2 ; 



3«i + |ci(u + u); 



«' = k + v^ 2 + ; 
\/ v 2 + ui, v' = v — k; 



v 7- " + ^1 

I (±u? - U2 + -wi + cie" + 2c 2 e 2u ) 
y/2u, v 



\v 1 + ciVit + vi; 



t/ = f («i + cie u ); 



-§W2 



^ 2 



cV 2 ^; 



2ciU, v' = - 
V3u ie v + 2ci 



■§«2 



3 7 ,2 



in 2 e 2 ^ 



3^2 



-2u 



the list 


of Sec- 


connected with 


m - 




m - 


-» 


m - 


()3.10|) 


m> - 


• (13.111) 


m - 


(j3.12[) 


m - 


• (j3.13|) 


m - 


• (|3.14|) 


m> - 


• (3.16} 


m - 


• (]3.17() 


<m - 


• 03.18[) 


m - 


• (|3.22|) 


m - 


• (j3.23|) 


m> - 


■ (13.24)) 
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v' = | v 2 - | v\ + 2V3c lVl e- v - \ e 2v (u 2 + c 2 ) - c\e~ 2v - ^ -» 

u' = §\Z2V(3«i + cie 1 " + 2uV), 

u' = 3 U2 -| t ,2_ Uie « + i e 2v ( Cl - 2u 2 ); P^Sjl 

/ _ 3 
2 



u 



u 



u' 



V2e v {ui + 2c 2 u + cxu 2 e v ), 



v' = | v 2 - I vf - lc lU ie v + 3c 2 d - |(2c 2 + ciue") 2 - ft^e 2 "; fl3U -» (E2SD 



' — 3 
2 

/ _ 3 



V2e v ( Ul + (c + u 2 )e v ), v' = \v 2 - \v\ - lu ie v + \{c - u 2 )e 2v ; -> (MI) 

\/2 (!*2 - ui«i) + \/2 (c 2 e" - c 3 e~ M - 3ciU X e -1 ') , 
v > = lv 2 -l(ul + vf) + c 2 e u + c 3 e~ u - 6cim e - " - 3c? e" 2,; ; ([3U) -» <|532j) 

ii' = I\/2(« 2 - UiUi) - 3\/2(cui + 2k Vl ) e~ v - 6ckV2 e~ 2v , 



2 



2 

u' = | U2 -|(^ + w 2 ) - 3(fe«i + 2cui) e _ " - 3(c 2 + A; 2 ) e" 2 "; JHJ -» (1333]) 



^ (ua - ui«i) - \/2 (c2 e~ u + 6cic 3 - 3ci«i e u+v + 3c 3 ui e - ") , -» (MI) 



r/ = | ^ 2 - f (n 2 + w 2 ) + c 2 e~ u - 3cj e 2{u+v) - 3d«i e u+v - 6c 3 d e~ v - 3c 2 , e~ 2v ; 
v! = -4 (u 2 - ui«i) + 3\/2 fci e 2 ^) + c 2 ui e u+w 



t/ = I W2 - f (u? + vj) + 3( Cl - c|) e 2 ^ - 3c 2Ul e u+v ; (33]) - (Ml 

u' = x/2(| (u 2 - uid) - | fcc 2 e" 21 ' + 3c 2 m e u+v + 2kc x c 2 e u ) 
+ dV2 (km + 2v 1 )e" v , k 2 = l, 

v l = l V2 -3( u l+ v j) - | c \ e ~2v _ 3c 2 e 2(u+v) _ 3^ & u + v 

+ 2dc 2 e u + d(«i + 2kvi)e- v ; JUJ -» <|53Sj) 



/ _ 3 
2 



V2 (u 2 - uid) + 3^2 (d e u+v + c 2 e"-" - c 3 e" 1 ') «i 



+ 6c 3 ^/2(c 2 e- u -d e u ), 

v '=3 V2 _3 {u 2 + v 2 } _ 3c 2 e 2(u+v) _ 3^ ^-u) 

+ 3ui(c 2 e" - " - d e n+,; ) - 6c 3 d e~ v - Qc x c 2 e 2v - 3c 2 3 e~ 2v ; -> (|5^7|l 

u = § V2 (u 2 - uid) + 3V2 e" 2 ^ - i V2 c 2 ( Ul - 4d) e" 1 ' 

- 2V2 dd e u + § y / 2d(2d e 2(u+t,) + 3«i e"^), 

/ 3 3^2, 2\ 15 2 2(u+v) 9 u+v , 5 u 

- ^e-^ + c^m -d)e" ,; ; P3]l -> (^Sjl 
«' = § y/2 {u 2 - «xd) + 3^2 ui (d + c 2 e 11 "") - §v^c§ e~ 2,; 

+ 2^2 c 3 (c 2 e- u - d e u ) - \/2 c 3 (ui + 2v 1 )e~ v , 

t/ = | U2 _3 (u 2 + y2) _ 3e 2, (ci e u + ^ e - U) 2 _ Cs(ui + 2wi)e -, 

+ 2c 3 (c 2 e- u - d e u ) -14 e~ 2v + 3m (c 2 e""" - d e M+,; ) ; (JO 



^" 2 " 3 ^" 2 2g3 - 3^2 (3, + 2c)V^7, 



v =3(d -2c 5 -c 2 ), g = u + v; (f3TTI) -> (l3^0|) 

«' = u, t/ = ui; D d3I5D 



u' = \u 2 + u 2 + v, v = u, (1331) -»■ dSSD 



this substitution is invertible, see below (|3.6[) — > (j3.3|) ; 

u=ui + it; 2 , v ' = v-k; (|5^ j> (|5^| l 
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v! = 4u 2 + ui (f^u + vi - 2ci) + l§u 2 + \c\u, 



v' 




-3&-f«iW« + «i; 


(|3T3j) - 


+ (|3.18[) 


u' 




u\, v' = v; 


(1^4)) 


— > l|o.o|l 


11 




Ui, v'=vr, 






v! 




|u3 + 2uu\ +V\, v' = u, 


(ED 


-> (ESI) 

\) — : — y 


this substitution has the quasi-local inverse substitution, see below (j3.6[) - 


(HI; 




i 

u 


= 


U2 + vvi, v = v — k; 


jap 


-> 


i 

a 


= 


3n3 — 6(uiVi + VU2) + 2uu\, v' = 3vi — 3v 2 + u; 


o - 


+ (|3.14[) 


i 

a 


= 


e~ v , v' = \v 2 - \v\ - ue 2v ; 


d53D - 


* (|3.25|) 


a' 


= 


ui + V3e v (u 2 + Ui^i), 






i 

V 


= 


3 3 2 2v 1 AT —v 1 -2u 


- 


* (|3.26|) 


u! 


= 


v, v = -§f 2 - v 2 + u; 


dill) 


-» (USD 


a' 


= 


v, v' = —\v2 — v 2 + W, W x = U, Wt = U2 + uv; 




(E3D 


v! 


= 


vi, v = -%v 3 - v{ + u; 




-» 


f 

U 


= 


1 1 / 12 
u + k, v = vi — 2 v ; 




-> (1321) 


l 

a 


— 


v — k, v = u\ — ^V2 — 2~{v + 6k ) + zkv; 




-» (EH) 


v! 


— 


3u! - |ui - f(u 2 + u 2 ), 






1 

V 


— 


Q Q / ,f-»\,Q Q 2 , Q 2/n, \ 

|u 3 - ju(u 2 + 2v 2 ) + |w 2 - lu{ + |^(2ni + vt) 










-| Wl (2n + 7;)-^(n 2 - W 2 ) 2 ; 


dSSD - 


* (|3.10|) 


1 

u 




-3ii 2 + 






v' 


= 


|«3 + |«l«2 - - 2u 2 + 2 Vl ) 2 + |(n 2 + t-x) 2 - |4; 


(1331) - 


* ()3.1l!) 


u! 




-3ui + -u, 






v' 


= 


| U2 + * UVl -^(u 2 - 2u x + 2v) 2 + \(u 2 + v) 2 - |n 4 ; 


M - 

^ V 


+ (13.121) 


u' 




-3(ui + 2m) + |ci5, g = u + v, 






1 

V 


= 


|u 3 + |ci(u 2 - v 2 ) + I85U2 - 9n? + 9vf 










+ 9(c l5 + 2g 2 ) Ul - | c 2 (m + ff 2 ) + 36g 2 Vl + 9/; 




■* (13.131) 


u 




3fi + u — 3f 2 , 






v' 




— 1^3 + i^2 + 9u^ 2 — Qu\v + 6(3v 2 — + 3-y 2 (2n — 3v 2 ): 


dSSD - 


+ (|3.14|) 


,,/ 




^1+^2 . ; : 4 

+C1VW + -U1 oVl, 






v> 




3 «3 + V4 3 c ti2+^3 ( 8 u U l+ 7; 2 9 (""1 + V 2 f (u\ + V 2 ) 2 
2 ^U + Vl 4 y/u + Vl 3 a/u + V\ S (u + Vl) !3 / 2 U + Vl 
til + V 2 , , \\ 1 

8/ l 6 ( n 2 + U3) ci(ui + -y 2 )J + 2(-y 3 + 2n 2 + civ 2 ) 
(u + V\Y^ 










+ if (« 2 - «f) - |(f2 — «i — ci?;i)\At + ^i - ^c 2 (u + 2ui); 


llo-ol) - 


-> do. loll 


v! 




\/m + — k, 






v' 




+ (B1 + ^ + 2tv __ _ _ „ 

V^ + fi (m + 2 


M - 


- (|3.19|) 


u' 




iV6ui + u 2 + v, v 1 = |(u 2 — f 2 ) + iV§uv\ — v 2 ; 


dSSD - 


* (|3.20|) 



20 



A.G. Meshkov and M.Ju. Balakhnev 



u = 3(ui + v\ — ku), 

v' = -|w 3 - §uu 2 + 9(u + k)v 2 - f u\ - 9vf 

9 
1 



U 



+ 9u{u + k) Ul - 9{u 2 + k 2 )vi - \u 2 {u 2 + 2k 2 ); <^ -> (EOT!) 



u' = -3u 2 + |ui + |e w (ci + 4-^Ui), 

v ' = |«3 - | (5M? - «2 + ui + cie n + 2c 2 e 2u ) 2 + |ci« 2 e u + fui^ + \v\ 

+ |n 2 (n 2 + 4ui) + §e u (ci(5u? + 2ui) - V^(y 2 + 2ui^ + u?)) + 18c 2 2 e iu 



ci(ci - 12v^«i) + 4c 2 (3« 2 + 2«i)) + 18c 2 e 3 "(ci - 2 x /5 2 ~u 1 ); ([321) -» 



u' = -3u 2 - §«2 - I (u 2 + w 2 ) - 2^/6^ (ui + v x )e u+v - c X (T lv 

_ C2e 2(-+-) +C3e 2(,-«) ) 
«' = |U4 - |«1«3 + I «3(2«1 + Ul) - |ti 2 (ti 2 + lif) - - v\f 

+ §w 2 (u 2 + 1uxv x - v\) + 4^6^ Vl (2c 3 e 3v - u - c x e u - v ) + c 

+ 3-^/6^ (i> 3 - u x u 2 + ^ 2 (2ui + ui) + 2u\vi{ui + v{j)e u+v 

+ 6t; 2 ( C2 e 2 ^ + c 3 e 2 ^ - c x e~ 2v ) + \c x {u x - ui)(«i - S^e" 2 ^ 

- c^e" 411 + 4c 2 c 3 e 4u - §ca(tti - ^) 2 e 2 ( M+ ^ + |c 3 (n 1 - 3^) 2 e 2 ^; S -> 
u > = 3u 2 - \v 2 - \ (u\ + v\) - 3ke~ v ( Ul + 2vi) - 3(ci + 2£; 2 )e- 2l \ 

t/ = |u 4 - |ui« 3 + |«3(ui " 2ui) + |u 2 (u 2 - u 2 ) - ^(u 2 - vl) 2 
+ |w 2 (n 2 - 2«it;i - v\) - 9cle~ Av + ISkcie'^ {u x + 2vt) 
+ |e" 2l, (4/c 2 ti 2 - Aciv 2 - 2k 2 (m + vi) 2 + ci(m + ui)(m + 3vi)) 

- ffce"" (u 3 + 2v 3 - Au 2 {u x + v x ) + v 2 (2vi - 3m) + u\(ui + 2«i)) , 

Cl = c 2 - A; 2 ; -> d£MD 



3« 2 - |«2 - |(n 2 + w 2 ) - 3c 2 e 2 ( M+ ^ - 3ci(3«i + 2 Vl )e u+v 



2 WZ 4 

+ c 2 e~ u - ke~ 2u , k = 3c 2 3 , 
v = f(t> 4 - u x u 3 ) + jv 3 (2ui + v x ) - \u 2 {u 2 + u 2 ) + |v 2 (u 2 + 2u\V\ - vl) 

_ ^( n f _ ^2)2 _ 9c 4 e 4( U+ ,) _ 18c 3 (3ui + 2vi )e 3 ^ + 6 C ? C2e " +2u 

+ 6kcje 2u + 2A;c 2 e- n - 21 ' + 6ka(ui - 2v x )e u ~ v + \c 2 eT u (u 2 + 2t; 2 + v\) 

- lde u+v (u 3 - 2v 3 + u 2 (9m + 4t7i) - u 2 (iti + 2«i) + 4n 2 (ni + Wl )) 

- \ke~ 2v {\v 2 +ui(4u 1 - ui) -3u 2 ) - fc 2 e- 4v + 6ci ^"(m + 2«i); -> (BHD 
u ' = 3n 2 - \v 2 - \ (u 2 + t; 2 ) + 3(ci - c|)e 2 ( u+u ) + 3c 2 ( Ul + 2 Vl )e u+v , 

9„, 9„, „, 1 9„, /„, n„. \ 1 9„, /„ ,2 „, \ 9 /„,2 2n2 



v = |u4 - i«i^3 + 3^3(^1 - 2ui) + |u 2 (wj - U 2 ) - fg(Ux - uf)' 

+ |w 2 (n 2 - 2u lVl - vl) - 9kc 2 2 e 4{u+v) + 18/tc 2 e 3(n+,;) (ni + 2v x ) 

+ |e 2(u+1,) (2c|n 2 + 4 Cl t; 2 + ci(«i + 3vi) 2 - t|(u? + 4ui«i + 5v 2 )) 

-|c 2 e ti+, '(u 3 + 2?; 3 +u 2 (n 1 +4v 1 )+?; 2 (27;i-ui)), fc = c| - 2cr, QMS) ^ QSM 
u > = 3u 2 - \v 2 - \ (u 2 + vl) - 3c 2 2 e 2 ^ + 3c 2 ( Ul + 2v x )e u+v 

+ 2c lC2 e u + Cl e- W («i + 2v x ) - \c\e~ 2 \ 
v = |u 4 - |iiiii3 + |w 3 (ui - 2«i) + |u 2 (wf - u 2 ) - ^(u 2 - f 2 ) 2 

+ |w 2 (n 2 - 2«i«i - u 2 ) - 9c^e 4(u+,;) + 18cle 3 ( u+,;) ( Ul + 2«i) 

+ 12cicie 3u+21 ' + |c|e 2(tt+1,) (2n 2 - u 2 - 4«ii^ - 5« 2 ) 

-4 C 2 C 2 e 2 " + C 2 e - 2 -(2n 2 -(n 1 + Wl ) 2 ) 
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- %c 2 e u+v (u 3 + 2v 3 + 1*2 («i + 4ui) + v 2 (2^i - m)) 

- Qcic 2 2 e 2u+V {ui + 2ui) - 3cic 2 e u (3ii2 + 2u 2 + 2t*iui + 3v\) 

+ fcie" 1 ' (« 3 + 2v 3 - 4-u 2 (ni + vi) + «2(2t7i - 3«i) + + 2«i)) ; ([221) -> (133BD 

u ' = 3n 2 - i«2 -|(«? + «?) - f c\e 2 ^ + | Cl ( Ul + 2^ 1 )e u+u 

+ \c x c 2 e u + c 2 e _,; (ui + 2v x ) - ^c 2 2 e~ 2v , 
v' = |u 4 - §«i« 3 + |u 3 (ui - 2i*i) + |« 2 (uf - u 2 ) - j^{u\ - v 2 ) 2 

+ lv 2 (u 2 - 2u lVl - v 2 ) - fief e 4 (« + ") + f c? e 3 ^') (t*i + 2^) 

+ f C f C2 e 3 «+ 2 " + | c 2 e 2(«+,)(l 8u2 + 16v2 _ 5u 2 _ Uuin _ 
27 2 2 2m . 3 3 u— 2v 9 2 u— 1>/ , o \ 

-f-c^e +fcic 2 e -|cic 2 e (u 1 + 2vi) 

- ^-de u+v (u 3 + 2v 3 + u 2 (m + iv^ + v 2 (2 Vl -ui)) 

- \ Cl c 2 e u {l7u 2 + 14w 2 + 12ui«i + 19u?) + £cie -3u (ui + 2vi) 
+ \cle- 2v {lQu 2 + 12v 2 - llu\ - 28nivi - I7v 2 ) - ^4 e ~* V 





+ |c 2 e v (u 3 + 2vs - 4u 2 (ui + v\) + ^2(2^1 - 3ui) + u\(u\ + 2v\j) 


; iD - 


>■ (13.381) 


u 


= 3— — U2 + 2gVl _|_ e 5v / Ul _ g 2 _ q Ui _ 3ui _ 3c (2 9 + c ) ) 5 = 

V «i - r 


w + 




1 

V 


= §U3 + 9(3 5 + c)u 2 + 9(2(7 + C ) U2 - 9 (V ^2 + 2gt>i) 2 + g + ^ 2 

i*i - g l 






- I8u 2 + 18 V«i - <? 2 («2 + 2(2 5 + c)(vi + </ 2 )) + 36(2 9 + c)(«i - 5 


2^3/2 






- 9 ff 2 ni - 18(0/ + c) 2 + 2</>i - 36c<?((<? + cf + g 2 ) - 45 5 4 ; 


IIO.UII 


> (|3.4()|) 


v! 


= v-2k, v' = u - |ui + |A: 3 t; 


TO - 




this substitution is invertible: 






u' 


= |«i + u - |/c 3 t, u' = u + 2A;; 


n - 


(EH) 


i 

u 


= — u' = yju + «i; 


- 


. jo inii 

» (|o.iyp 


v! 


= Cl\ / 2e" , ', v' = -vi + \ue v ; 


(|3.10|) - 


> (|3.24j) 


1 

u 


= ue , v = v\ — c\ue — 2c 2 ; 


(|3.10|) - 


> (|3.29[) 


u 


= (u + \/— c)e 1 ', u' = vi + (v— c — u)e 1 '; 


(|3.10j) - 


>■ (13.301) 


v! 


= iti + 2A;e-^, v' = v x - 2ce~ v ; 


(|3.10|) - 


> (|3.33[) 


u 


= Ul + ce u+v , v' = Vl -ke u+v , 






where c and k are roots of z 2 — 2c 2 z + 2c\ = 0; 


(|3.10|) - 


► (|3.35|) 


u 


= ui - fci e" v , «' = ui - 2c 2 + §ci e""; 


(|3.10|) - 


> (|3.36[) 


u' 


= ui + 2ci - | c 2 e-", i/ = «i - ci e M+,; + | c 2 e""; 


(|3.10|) - 


> (|3.38|) 


v! 


= 2^ ui — (u + v) 2 , v = 2(u + v + c); 


(|3.10|) - 


(|3.40|) 


v! 


= -2u, u' = |n-2(wi + *; 2 ); 


(]3.12j) - 


> (|3.14|) 


this substitution is invertible: 






v! 


= |(u -i*i) + fi* 2 , v ' = -\u\ 


(|3.14j) - 


> (]3.12|) 


u' 


= \{u\— uv), v' = — u); 


d3.14j) - 


> (|3.1()|) 


u' 


= 3(it + u) 2 - 3ui + ^c 2 , v' = \c\-u - v; 


(|3.14[) - 


► (|3.13[) 


u' 


= u, v' = V\\ 


(|3.14|) - 


> (|3.15|) 


u! 


A 3 2 / I 2 / 


(|3.14|) - 


► (|3.18|) 
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u, v 
civ2u, 



l v l + 77? ° 



V 

1 2v 



2 



- Ul e v - i e 2t ' (2n 2 + Cl ) , 
uie v + | e 2v (2u 2 + ci) , 



= — |«i + § ue"; 



3„D 

2 

9 r(u + U 



2 „u\ 



e v (u\ + 2c2U + c\u e 

2 9 



■iui + i(ci + l)ue v + C2; 



this substitution is invertible: 



1)3.14} - 


■* (|3.23|) 


(13.141) - 


■+ (13.241) 


(13.141) - 


-> (13 281) 


(|3.14[) - 


+ (|3.28[) 


(|3.14j) - 


+ (|3.29[) 


(|3.18[) - 


+ (13.131) 



9 u i 



V 



([TOD -> (I3~T8|) 



The graph of the substitutions is very cumbersome, therefore we show the most interesting 
subgraph only. 




Fig. 1. A subgraph of the differential substitutions. 



Comments. (1) System (|3.ip coincides with (|1.2p and (|3.3p coincides with (jl.3p . 

(2) It was simpler to obtain some systems from (]3.6[) than (I3.3P or vice versa. These systems 
are connected by the second order invertible substitution, see (|3.6p — > (|3.3p and f|3.3[) — ► (|3.6p . 
Hence, each system obtained from (13. 6|) can be obtained from (13. 3|) and vice versa. (3) Fifteen 
systems (^9]) - ([3~T5l) . (I3~T8|) . ([j\22]) . (f3T33]) ~ (f3736| ) . (1335]) and (foT40D can be obtained from both 
(|3.ip and (j3.6|) by the presented differential substitutions. Twelve systems (|3.2|) . (j3. 16|) . (j3. 17|) . 
d£Z3H , (pT241) . (I3~271) - (^30l ). (pT32]) . (pQ7j> and (pT39]) can be obtained from (J3HJ). The remaining 
eleven systems P^]) - (^3]) . (13T7D . (|5^j> . (1319]) (I3T2TT ) . ([3T25D . (IqT26D and (13T3TD can be obtained 
from (US} or (jOjl . 

Remark 7. As the systems (13. ip and (I3.3D have the Lax representations, then all systems from 
the list have the Lax representations in a generalized meaning (see Section [5]) . 

Remark 8. Some of the presented substitutions are superpositions of lower order substitutions, 
other substitutions are irreducible. 

Remark 9. System (13. ip admits the first and second order substitutions and does not admit the 
third and fourth order substitutions. Probably it does not admit any higher order substitutions, 
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either. Systems fj3.4j) and f|3.6|) admit the substitutions from the first till fourth orders. We 
do not present higher order substitutions for (|3.4|) because simpler substitutions exist for (|3.6|) . 
Fifth and higher order substitutions for systems (|3.4|) and (|3. 6|) have not been computed because 
the computations are extremely cumbersome. 

Remark 10. There are some additional differential substitutions under the constraints for 
constants in the systems. For example, there are substitutions (|3.ip — » (|3.3ip for c\ = and 
(|3.6|) — > (|3.32p for C3 = or (|3.6p — > (|3.39p for C2 = and so on. These substitutions are not so 
important and we do not present them here. 

Remark 11. Unexpectedly, the well known systems (I3.ip and (|3.3p are implicitly connected as 
it is shown in Fig. 1. 



5 Examples of zero curvature representations 

The 1ST method for nonlinear equations with two independent variables is based on investigation 
of a linear overdetermined system 

L(u,X,d x )^ = 0, ip t = A(u, A, d x )ip, (5.1) 

where L and A are ordinary linear operators, u is a smooth (vector) function satisfying a non- 
linear partial differential equation E(u) = and A is a parameter. The operators L, A may be 
both scalar and matrix. The operators L, A and E may be also pseudodifferential or integro- 
differential. The compatibility condition of system (|5.ip reads 

% +LA )n =°- 

Ot J I Lip=0 

There are two ways to satisfy this condition. The first operator condition was introduced by 
P.D. Lax [37]: 

FIT r)T 

— + LA = AL, or ^ = [4,4 (5.2) 



The second more general operator condition was introduced in [38], see also [39 : 

dL , s 

— + LA = BL. (5.3) 

If an equation E(u) = is equivalent to equation (|5.2p . then (|5.2p is said to be the Lax 
representation of the equation E(u) = 0. The pair of operators (L, A) is said to be the (L,A)- 
pair or the Lax pair. 

If an equation E(u) = is equivalent to equation (]5.3[) . then (15. 3p is said to be the (L, A, B) 
representation of the equation E(u) = or the triad representation. 

In all cases, operator L must essentially depend on the parameter A. This parameter cannot 
be removed by a gauge transformation L — > f Lf with some smooth function /, in particular. 

If system (|5.ip is differential, then the standard substitution tp = $1, ip x = ^2 and so on, 
provides the following first order system: 

^ = UV, % = VV, (5.4) 

where U and V are square matrices depending on u and A. The compatibility condition of linear 
system (|5.4p reads 



u t - V x + [U, V] = 



(5.5) 
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if E(u) = 0. Usually a stronger condition is required: (|5,5p is valid iff E{u) = 0. In this 
case equation (|5.5|) is said to be the zero curvature representation. For an evolutionary system 
ut = K(u,u x , . . . ,u n ), u = {u a } the matrix U usually depend on u only, but it may depend 
on u, u x , u xx , and so on. Let us consider the general case. 

If some smooth functions F = F(u, u x , ■ ■ ■ , u r ) and <I> = &(u, u x , . . . , u p ) satisfy the condi- 
tion 

(d t F + $)\ = 0, (5.6) 



then one obtains 



dF OF OF 

ft d t F\ = T—d t uf\ = T—dtu?\ = T—DI.K 



I u t =K \u t =K OU i \ut=K O u i I u t =K "-j 

where the summation over i = 0, . . . , r and a = 1, . . . , m is implied. This implies 
dF 

$ = — - — DlK a 
duf x 



according to (|5.6p . Using this result one obtains the following identity: 

d t F + $^-^D x (u?-K a ), \fu (5.7) 



for any F and $ satisfying (|5,6|) . 

Let us apply identity (|5.7p to equation (15. 5p . If the matrix U depends on u only then 



It is obvious now that equation (|5,5p is equivalent to Ut = K iff the matrices dU/du a , a = 
l,...,m are linearly independent. Suppose now that the matrix U depends on u and u x then 
one obtains 



If the matrices 8U/du a , dU/du x , a, f3 = 1, . . . ,m are linearly independent, then equation (|5.5|) is 
equivalent to Ut = K again. Otherwise, equation (|5.5p would be equivalent to some differential 
consequence of the system Ut = K that is a more general system than the original one. In this 
case we call equation (|5.5p the generalized zero curvature representation. 

It is well known that equations (|5.4|) and (15.5f) are covariant under the following transforma- 
tion 

y = s~ 1 y, u = s-\us - s x ), v = s~ 1 (ys-St), (5.10) 

where S is any non-degenerate matrix. This transformation is called a gauge one. Any gauge 
transformation is invertible and preserves compatibility of system ()5.4p . 

Two (L,A)-paiis were proposed for system (jl.2p in [3]. One of these (L,74)-pairs coincides 
with the (L, A)-pair that was presented in [4j. The L-operator of the common (L, ^4)-pair takes 
the form L = (d x + f){d x — g), where 

f = UuV2-2v), g = UuV2 + 2v). 
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The temporal Lax equation reads ipt = Aip, where A is a fractional degree of L. The spatial 
Lax equation Lip = X 2 ip can be transformed into the system (d 2 — g)ip = Xcp, (d 2 + f)ip = Xtp 
and then into the normal form (15.41). were 



U 



1 

(n v / 2-2u)/6 



A 







V 




A 

1 

uV2 + 2v)/6 Oj 

/1 + /2 -(«iv^ + «i)/6 Xv/3 



-2A 




\ 



V 





Xv/3 



-2A 




Here f\ and f 2 take the following form: 



fi = ±(3v 2 - 2u 2 + 2v 2 



2A 2 , 



i Ul V2- Vl )/6 { Ul V2- Vl )/3 
h-h ( Ul V2- Vl )/6j 



h = TF(3^2 +uv). 



(5.11) 



Matrices (|5.1ip realize the zero curvature representation of system (|1.2p . 

System (I1.3P also has two Lax representations (see [3]). Using the simpler L-operator, we 
have found, similarly to the previous case, the following matrices that realize the zero curvature 
representation of system (11 .3[> : 



U 







/ 10 

-v/3-X 1 

1 

\ u/9 X-v/3 0/ 



V 



( ui/6 
h - Xv/3 
ui/9 
\uv/27 + u 2 /9 



2X-v/3 
-vi/6 
-2u/9 
-ui/9 





v/3 
vi/6 
h + Xv/3 



\ 







-2A - u/3 
-«i/6 / 



(5.12) 



where 



h = \v 2 + \{v 2 - 2u) - 2X 2 . 

Let us consider an admissible differential substitution u = f(ui,Vj), v = g(Hi,Vj) of sys- 
tem (jl.ip . Substituting u and v in the matrices U(u,v) and V(ui,Vj) one obtains 

U(ui,Vj) = U(f,g), V(ui,Vj) = VUtf.l )',!,). 

As U depends on derivatives of u or v, then one has a generalized zero curvature representation. 

To obtain an ordinary zero curvature representation one can try to remove higher order 
derivatives from the matrix U using the gauge transformation (|5.10p . But this is not always 
possible (see example B below). 

A. Performing the substitution (|3.ip — > (|3.9p into the matrix U from (|5.1ip one obtains 



U 



I 1 

(u 2 -v 2 )/2 + h 1 



\ A 




A 




0\ 



1 



-{u 2 + v 2 )/2 + h 2 0/ 
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where 



hi = \{ui - vi) 2 , h 2 = \{ui + vi) 2 . 

One can easily verify that matrices A = U U2 and B = U V2 are commutative, hence the system 
S x = (Au xx + Bv xx )S has the following solution S = exp(Au x + Bv x ). The matrix U\ = U 
evaluated according to (15.10P takes the following form 



(u x -\ 



V 2A 



v x - u x 






2A 

u x + v x j 



Now another gauge transformation is possible with the following diagonal matrix: 



S\ = exp diag (Ui)dx^j , 



where diag([/i) is the main diagonal of Ui. This gauge transformation provides the following 
matrix 



A corresponding V^-matrix can be obtained by solving equation (15. 5p directly or by the previous 
twofold gauge transformation. This matrix takes the following form 



f 


e v-u 








\ 








Xe~ v 


















V 













) 



( 



Vo 





-2X 2 e u ~ v 
Xe u (vi — u\) 

^e- v (2v 2 -3u 2 + v 2 ) 



where 



h 



VI 



-u 2 



e v ~ u fi 


-2\e v 
-\e- u { Ul + Vl ) 



\e- u {ui + vi) 



-2Xe v \ 



'{2v 2 -3ul + vf) Xe u ( Ul - Vl ) 





-2\ 2 e~ u ~ v 



e u+v h 





\v 2 + u\v\ + \{u\ + vf), h = u 2 - \v 2 - U\V X + \{u\ + vf). 



The matrices U 2 and V 2 realize the zero curvature representation of system (|3.9p . 
B. Substitution (13. 1|) — > (I3.17|) reduces matrix U from (15. lip to the following form 



U 



( 





1 





OA 


(k- 


v)/3 + R 





A 
















1 




A 


o ( 


k-v)/3- R 


0/ 



where R = ^2{u\ + v 2 )/6. It is obvious that one cannot remove u\ from U by a gauge trans- 
formation. It is clear from the structure of the matrix U that 

U t -V x + [U, V] = A {kvt + \u 2 - \ V?> - v t ) 

+ BD X (u 3 - f ( 2 ^ 1 + ^) + 3 VV2 + 3 V 2 + 2 y 3 _ k{ 2 v 2 + j _ 
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where A and B are some linearly independent matrices. Thus, this zero curvature representation 
for system (|3.17p is generalized. Of course, one may introduce here the new variable u' = u\ to 
obtain an ordinary zero curvature representation. But we do not know if it is always possible. 
C. Performing the substitution (|3.3p — > (|3.6p into the matrix U from (|5.12p one obtains 



U 



( 1 0\ 

-n/3-A 1 

1 

\u 2 /6 + (n 2 + v)/9 A-n/3 0/ 



The first gauge transformation is performed using 5*1 = exp(ui(dll /du2)): 



Si 



(I 0\ 

10 

10 

\m/6 001/ 



The transformed [/-matrix is 



( 1 0\ 

-u/3-X 1 

m/6 1 

\{u 2 + v)/9 -m/6 A-n/3 0/ 



The second gauge transformation is performed using S 2 = exp(n {dU\ / dui)): 



( 1 













1 








n/6 





1 





V o 


-u/Q 





V 



The result of the twofold gauge transformation is 



U 2 



Vo 



( 1 0\ 

-n/6 — A 1 

-u/3 1 

\u 2 /36 + v/9 A - n/6 0/ 



/ 


-ui/6 


2A 







-2\ 




h - An/3 


-ui/6 




u/3 





(uui 


- ui)/18 - Am/3 


-n 2 /3-n 2 /6- 


- 2v/9 


m/6 


-2A 


\ 


h 


(^i — UU\)/1S — 


Xui/3 


h + An/3 


ui/6/ 



where 



\u 2 



1 „.2 



IS 



u 



\v-2\ 2 . 



/ 4 = j|(w!/2 -v 2 + u\) + jggn 3 



27 



uv 



InA 2 . 



Matrices f/2 and V 2 realize the zero curvature representation of system (|3.6p . 



6 Conclusion 



The examples in Section [S] illustrate the fact that some systems possess ordinary zero curva- 
ture representation while others possess generalized zero curvature representation. All these 
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representations are obtained from the Drinfeld-Sokolov L, A operators by using corresponding 
differential substitutions listed in Section UJ Matrices U and V that realize all zero curvature 
representations have the size 4x4. Thus, the two-field evolutionary systems presented above 
are integrable in principle by the inverse spectral transform method. But the fact is that the 
inverse scattering problem for differential equations with order more than two is extremely 
difficult. That is why other methods for solution of equations may be useful [40]. They may be 
Backlund transformations [41 j . Darboux transformations |42U43| . Hirota method [H] or numeric 
simulating (see [45], for example). 
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